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Beschreibung 

Die Erfindung betrifft dn Verfahren fiir cioen Coniputcr- 
tomograpben (CT) zur Nachverarbeitung cines Schnitlbil- 
dcs sowie cincn nach dnem solchen Vcifahren arbcitenden 
Computertomograpben. 

Durcb djfferentielle Abweichungen dcr dnzelnen Mefisy- 
stemkanile eines CompuleriomographeD vom jeweiligen 
Kalibrierzustand cntstehen bei CT-GerStcn der 3. Genera- 
don in den Schnittbildern Ringartefakte, die uni das Dreh- 
zentru m des GerStes zentrieit sind. Dieses Problem ist in 
IEEE Proceedings, Vol. 134, Pt A, February 1987, Seiten 
126-135, beschriebcn. 

Dcr Erfindung liegt die Aufgabe zugrunde, dn Verfahren 
und dnen ComputertonnogFaphen der eingangs genannten 
Art so auszubilden, dafi eine Korrektur dieser Ringartefakte 
in den rekonstniimen Schnittbildern in Form dner Bild- 
nachveraibeitung erfolgt. 

Diese Aufgabe ist erfindungsgemafi gelQst durch die 
Merkmale der PatentansprQche 1 und 3. Bei der Erfindung 
erfolgt die Detektion der Ringartefakte durch Mcdianfilte- 
rung in radialer Richtung. Ringartefakte werden so vollstSn- 
dig konigiert und nicht lediglich in ihrer Amplitude rcdu- 
ziert und verschmiert Im Unterschied zu bekannten Verfah- 
ren werden femerhin stfirende Knochenstrukturen im Bild 
vor der Medianfilterung durch eine einfache Schwellenope- 
ration eliminiert 

Die Erfindung ist nachfolgend anhand der Zeichnung nS- 
her eriautert Es zdgen: 

Fig. 1 die wesentlichen Ibile dnes Computertomogra- 
phen der 3 . Generation, 

Fig. 2 eine Darstellung der Rechenschritte ftir den Com- 
putertomographen gemSB Fig. 1, 

Fig. 3 dne schematische Darstellung zur Filterung fiir das 
erfindungsgemSBe Verfahren, und 

Fig. 4 eine schematische Darstellung zur Wnkdvcr- 
schldfting fiir das erfindungsgemSBe Verfahren, 

Der in der Fig. 1 dargestellte Computertomograph weist 
eine McBdnhdl aus eincr RGntgenstrahlenquelle 1, die ein 
facherfbrmiges ROntgenstrahlenbUndel 2 aussendet, und ei- 
ncn Detektor 3 auf, welcher aus dner Rdhe von Einzelde- 
lekioren, z. B, aus 512 Einzeldetektoren. bestehL DerFokus 
ist mil 11 bezeichnet Der zu untersuchende Patient 4 liegt 
auf einer Patientenliege 5. Zur Abtastung des Patieaten 4 
wird die MeBeinheit 1, 3 um ein MeBfeld 9, in dem der Pa- 
tient 4 liegt, um 360^ gcdrcht Die Drchachsc ist.mit 10 be- 
zeichnet. Dabd wird die Rontgenstrahlenquelle 1, die von 
einem Rdntgengenerator 6 gespeist wird, gepulst oder mit 
Dauerstrahlung bctrieben. Bd vorbestimmten Winkelposi- 
rionen der Mefidnheit 1, 3 werden S5tze von Daten'erzeugt, 
die vom Detektor 3 einem Rechner 7 zugefiihrt werden, wel- 
cher aus den erzeugten DatensStzen die SchwSchungskoeffi- 
zienten vorbestimmter Bildpunkte berechnet und auf einem 
Sichtgerat 8 bildlich wiedergibt Auf dem Sichtger^ 8 er- 
schdnt demgem&fi ein Bild der durchstrahiten Schicbt des 
I^tienten. 

In der Fig. 2 siellt der Krds 16 schematisch das Eingangs- 
bild dar. In den Verfahrensschritten 17 bis 25 erfolgen nach- 
einander folgende Bildverarbeimngsschritte: 
Schritt 17 Abschneiden von Knochen und Luft im Bild; 
Schritt 18 Medianfilter mit Raster ARl flber ganzes Bild 
usw.; 

Schritt 19 Schwellwertbewertung; 

Schritt 20 Differenzbilderzeugung; 

Schritt 21 Zusatzb'ches Medianfilter b iiber mediangefilter- 

tes Bild mit Raster, 

AR2 im zentnmisnahen Berdch in einem Quadrat der Sei- 
teniangelAU; 
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Schritt 22 Schwellwertbewertung; 

Schritt 23 Differenzbilderzeugung; 

Schritt 24 Winkelverschleifting des Differenzbildes; 

Schritt 25 Subiraktion; 

. 5 Eingangsbild - winkelverschliffenes Differenzbild 
Der Krds 26 stellt dann das Ergebnisbild dar 
In dnem CT-Bild IBILD liegen Np • Np (z. B. 512 • 512) 
Pixelweite des Schnittbildes vor, in der Regel als 12 Bit In- 
tegerwerte im Wertebereich 0... 4095 ("Hounsfield 

10 Scala"), die ftir die weitere Bildbeazbeiuing auf die Integer- 
Werteskala-1024 bis +3071 umgesetzt werden. 

Im ersten Verarbdtungsschritt 17 werden Knochen- und 
Luftanteile im Bild eliminiert, um die sp^tere Ringdeiektion 
nicht zu stdren. Zu diesem Zweck werden alle CP-Werte 

15 groBer als dne bestimmte Scbwelle SWO (z,B. SWO = 
300), die von Knochenstrukturra stammen k5nnten, gleich 
SWO gesetzL AUe CT-Werte kleiner als eine bestimmte 
Scbwelle SWTJ, die von Lufl bzw. Lufteinschliissen stam- 
men kOnnten, werden gldch SWU (z. B. SWU = -300) ge- 

20 setzL Die Schwellen SWO und SWU sind sowohl vom ver- 
wendeten CT-Faltungskem, also der Maximalscharfe der 
CT-Rekonstruktion, als aucb vom maximalen SchwS- 
chungswert des daigestellten Objektes abhSngig und werden 
z. B. empirisch ermittelt. 

25 Dieser Verfahrensschritt 17 ergibt eine neue, mit RIBILD 
bezeichnete Bildmatrix, die ebenfalls Np • Np Pixdwerte 
hat. 

In den Verarbeitungsschritten 18 bis 20 wird ein Korrek- 
tuibild erzeugt, in dem die Ringartefakte isoliert sind. 

30 Zu diesem Zwecke wird in der Bildmatrix RIBILD zu- 
nlk;hst dne erste Medianfilterung voigenommen entlang ei- 
ner Vielzahl von durch das Drehzenttum des Computerto- 
mographen verlaufenden Geraden (radiale linien). Diese 
Geraden decken die Bildmatrix derartig ab, daB jedes Pixel 

35 auf dner solchen Geraden liegt Das erste radiale Medianfil- 
ter hat 2Mi + 1 Statzstellen (z.B. Mi = 2) im Abstand 
aai • aRi ist etwa halb so groS wie die halbe Linienbreite ei- 
nes Ringanefakts und ist daher sowohl vom gcwfihlten Aus- 
schnitt des MeBfddes (Zoomfaktor -Zoom") als auch von 

40 der Scharfe des verwendeten Fallungskems abhSngig und 
bei kleinem Bildausschnitt gr&Ber als der Pixelabstaod in ra- 
dialer Richtung. FUr jedes Pixel P im Bild RIBILD werden 
die im radialen Abstand -MiaRi, -(Mi+l)aRi, . . 
(Mi-l)aRi, MiBri liegenden Pixd z.B. durch nearest 

45 neighbour Interpolation besiimmt Der Pixelwert P des me- 
diangefilterten Bildes MEDIBILD ergibt sich aus dem Me^ 
dianwert Med der Pixelwerte der 2Mi + 1 Stiitzsteilen. Der 
Medianwert Med von 2Mi + 1 Werten ist deijenige Wert, fllr 
den Ml Werle kleiner oder gleich groB sind und die Ubrigen 

50 Ml Werte groBer oder gleich groB. Der Medianwert der 5 
Werte 1, 4, 3, 8, 17 z. B. ist Med = 4. Die erste Medianfilte- 
rung ist schematisch in Fig. 3 dargestellt. 

Nach der ersten Medianfilterung wird fUr alle Np • Np Pi- 
xel das Differenzbild DIFFl = RIBILD - MEDIBILD ge- 

55 bildet DIFFl soilte nur noch Ringartefakte enthalten. Um 
Reste sl5render Knochenkanten etc. zu eliminieren, wird im 
Bild DIFFl eine Schwellenoperation mit der Anefakt- 
schwelle voigenommen: alle Pixdwerte in DIFFl, die 
betragsm^ig grSBer als San sind, werden als StOiungen 

60 identifiziert, die Mschlicherweise als Ring erkannt wurden. 
Sie werden, wenn sie +Ssn Ubersteigen, gleich •^Sm und, 
wenn sie tmterschrei ten , -Son gesetzt San h^gt ab vom 
gewShlten CT-Faltungskem, d. h. von der Maximalscharfe 
der CT-Bildrekonstruktion, und wird z. B. empirisch ermit- 

65 telL Das mit der Artefaktschwelle San behandeltc Di fferenz- 
bild DIFFl hei6tR2BILD. 

Um eventuell vorhandene Reste von Ringartefakten mil 
anderem Rast^ zu eliminieren, wird fUr alle Pixdwerte 
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Np • Np das Differenzbild R3BILD aus dem maskicrten 
Ausgangsbild RIBILD und dem Bild R2BILD gebildeL Bei 
idealer Wiikung des ersten Mcdianfilters soUte R3BflJ) 
keinc Ringartefakte mehr enlhalten. In R3BEJ) (oder einem 
Ausschnitt davon) wird nun im Vcrfahrensschritt 21 eine 5 
zweite Medianfilterung voigenonuncn, wiedcnim enilang 
einer Vielzahl von durch das Drehzentrum (Drehachse 10) 
des Computertomographen verlaufendeo Gcraden (radiale 
Linien). Diese Geradcn decken die Bildmatrix derartig ab, 
dafi jedes Pixel des gewShlten Ausschnittes auf einer sol- 10 
Chen Gcraden liegt. Das zweite radiale Medianfilter hat 2Nt 
+ 1 ScUtzsteUen (z. B. M2 = 1) im Abstand ai^ (z. B. am = 
aRi/2). aR2 ist sowohl vom gewfihlten Ausschnitt des MeB- 
feldcs (Zoomfaktor "Zoom") als auch von der SchSrfe des 
verwendeten Faltungskems abhSngig und bei kleincm Bild- 15 
ausschnitt grbfier als-der Pixelabsiand in radialer Richtung. 
Fiir jedes Pixel P im gewUnschten Ausschnitt des Slides 
R3BILD werden die im radialen Abstand -M2aR2» -0^2 
+l)aR2, . . (M2-l)aR2, M2aR2 liegenden Pixel z. B. duich 
nearest neighbour Interpolation besdmnit Der Pixelwert P 20 
des mediangefilterten Bildes MED2BILD eigibt sich aus 
dem Medianwert Med der Pixelwerte der 2M2 +1 Stutzstel- 
len. 

Nach der zwdten Medianfilterung wild fUr alle Np • Np 
KxeldasDifrerenzbildDIEF2«RlBILD-MED2BILDge- 25 
bildet Um Reste stdrender Knochenkanton etc. zu eliminie- 
ren, wird im Bild DIPF2 wiederum eine SchweUeooperation 
mit der Artcfaklschwelle San vorgenommcn: alle Pixelwerte 
in DIFF2, die bctragsmiBig gr6Ba: als sind, werden als 
StSrungen identifiziert, die fdlschlicherweise als Ring er- 30 
kannt wurden. Sie werden nach der bereits beschriebenen 
Voigchenswdse gleich +S«n oder -S^it gesetzt. wird 
z. B, empirisch ermitlelt und hSngi wie erwShnt ab vom ge- 
wahllen CT-Fallungskem, d. h. von der MaximalschSrfe der 
CT-Bildrekonstruktion. D as mi t der Anefaktschwelle Stn 35 
behandeltc DifFerenzbiW DIFF2 heifit R4BILD. 

Im Verarbeitungsschritt 24 werden stfirende Rauschstruk- 
turen in R4BILD eliminiert. R4BILD enthait nicht nur Rin- 
gartefakte, sondem durch das unvenneidlichePixeliauschen 
in CT-Bildem auch Rauschstrukturen. Zur UnterdrQckung 40 
dieser Rauschstrukturen macht man sich zunutze, daB Rin- 
gartefakte auf Kreisbdgen einer gewissen Mindestlange 
mehr oder wehiger konstant sdn mlissen, um flberhaupi als 
Ringe erkannt zu werden. Der drittc Veraifoeitungsschritt ist 
also eine TiefpaBfilterung der Pixelwerte des Bildes 45 
R4BIID entlang von Kreisbogen um das Drehzentrum; da- 
durch bleiben die Ringartefakte erhallen, Rauschstrukturen 
dagegen werden weggeglStteL Die "nefpaBfilterung kann 
z. B. als einfache Mttelwertbildung realisiert werden: fUr 
jedes Pixd P im Bild R4BILD mit den Kooidinaten xp, yp 50 
sucht man auf dem durch P verlaufe nden Kreisbogen um das 
Drehzentrum mit Radius ^x^pfy^p die 2Nw + 1 Pixel mit 
dem Winkelabstand -Nw SPHI, -<Nw+l) • SPHI, . . . 
(Nw-1) • SPm, Nw • SPHI von R SPHI ist dn ftei wShlba- 
res Winkelinkrement Der Wert des Kxel P wird erselzl 55 
durch den Mittelwert der 2Nw 1 Pixel. Die Pixel werden 
dabei je nach Abstand vom Drehzentrum in Pixel im Innen- 
bereich und Pixel im Aufienberdch eingeteilt, um unter- 
scbiedliche ^nkelinkremente SPHI fUr die Mitteilung im 
Innenbereich und im AuBenbereich realisieiren zu kdnnen. 60 
Durch die Winkelverscbleifiing eatsteht aus R4BILD das 
Bild WBILD. Die >^^nkdverschleifung ist schematiscb in 
Fig. 4 dargestellL 

Subn^ktion von IBILD und WBUD im Schriu 25 
schlieBlich ergibt das Ergebnisbild, in dem die Ringartefakte 65 
besdtigtsind 

Die Fig. 3 zeigt dne schemadsche Darstellung der ersten 
Medianfilterung gemaB dem iSchritt 18. 
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Die Fig. 4 zeigt dne schematische Darstellung der Win- 
kelverschlcifung gemaB dem Schriu 24. 

Erfindungswesentlich ist, daB ein Schnittbild des gesam- 
ten MeBfeldes oder nur eines Ausschnittes daraus in der 
Weise nachverarbeitet wird, daB die Pixdwerte des Schnitt- 
bildes einer oder mebroren Medianfilterungen und einer 
Mictelung unterzogen werden, wobd die Medianfilterungen 
und die Mittelung entlang einer Anzahl von AusfUhrungs- 
richtimgen erfolgen. 

Patentanspriiche 

1. Verfahren fiir einen Computertomographen mit ei- 
nem Universalrechner (7) zur Nachveraibeitung dnes 
rekonstruicrten Schnittbildes einer Schicht eines Un- 
tersuchungsobjektes (4) zur Beseidgung von Ringarte- 
fiakten, aufweisend folgende Verfahrensschritte: 

a) Maskierung von Knochen- und Luftanteilen 
im Bfld IBILD mit Np • Np Pixeln, wobei alle CT- 
Werte in . der Pixelmatrix groBer als eine be- 
sdnunte Scbwelle SWO gleich SWO gesetzt wer- 
den, alle CT-Werte kleiner als eine besdmmie 
Schwdle SWU gldch SWU, so dafi sich eine 
neue, mit RIBILD bezeichnete Bildmatrix etgibt, 
die ebenfalls Np • Np Pixelwerte hat. 

b) Durchfuhning dner Medianfilteroperadon in 
RIBILD, die vorgenommen wird entlang einer 
Vielzahl von durch das DrehzenUum (10) des 
Cbmputertomogr^hen verlaufenden Geraden 
(radiale Linien), wdche die Bildmatrix derardg 
abdecken. daB jedes Pixel auf einer solchen Gera- 
den liegt, wobei das Medianfilter 2Mi + 1 Stiitz- 
stellen im Abstand aRi hat und der Pixelwert P des 
mediangefilterten Bildes MBDlBiki sich ergibt 
aus dem Medianwert Med der Pixelwerte der 2Mi 
+ 1 Statzstellen von RIBILD: 

c) BildungdesDifferenzbildesblFFl ss RIBILD 
- MEDIBILD fiir alle Np - Np Pixelwerte und 
DurchfUhrung einer Schwellenoperation mit der 
ArtefaktschweUe S^rt im Bild DIFFl, wobd alle 
Pixelwerte in DIFFl, die bctragsmaBig groBer als 
Sart sind, gleich +Sart oder -San gesetzt werden, so 
daB das Bild R2BILD entsteht. 

d) Duichfijhrung einer TiefpaBfilterung der Pi- 
xelwerte des Bildes R2BTT.D witlang von Kreis- 
bogen um das Drehzentrum, die als einfache Mit- 
telwertbildung realisiert werden kann, wobei man 
fiir jedes Pixel P im Bild R2BILD mit den Koordi- 
naten xp, yp auf dem durch P verlaufende ai Kreis- 
bogen um das Drehzentrum mit Radius ^x^p+y^p 
die 2Nw + 1 Pixel mit dem Winkelabstand 
-Nw • SPHI, -(Nw+1) • SPHI, . . . (Nw-1). SPHI, 
Nw * SPHI von P sucht und den Wert des Pixels P 
ersetzt durch den Mittelwert der 2Nw -t- 1 Pixel, so 
dafi das Bild WBILD entsteht 

e) Subn^on von IBILD und WBILD, wobei 
das Ergebnisbild entsteht, in dem die Ringarte- 
fakte beseitigt sind. 

2. Verfahren fiir einen Computertomographen mit ei- 
nem Universalrechner (7) zur Nacbveraibeitung eines 
rekonstruicrten Schnittbildes einer Schicht eines Un- 
tersuchungsobjektes zur Beseitigung von Ringartefak- 
ten, aufwdsend folgende Verfahrensschritte: 

a) Maskierung von Knochen- und Luftanteilen 
im Bild IBUD mit Np • Np Pixeki, wobei alle CT- 
Werte in der Pixelmatrix grbBer als dne be- 
stimmle Scbwelle SWO gleich SWO gesetzt wer- 
den, alle CT-Werte kleiner als eine bestimmte 
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ScbweUe SWU gleich SWU, so daB deb eine 
ncue, mil RlBILD bezeichncie Bildmatrix eigibi, 
die cbcnfalls Np • Np Pixelwcrie hai. 
b) Durchftih];yng einer ersten Medianfilteropera- 
tion in RIBUD, die vorgenommen wird entlang 5 
einer Vielzahl von durch das Drehzentnim (10) 
des Computertomographeo verlaufcnden Geraden 
(radialen linien), welcbe die Bildmatrix derartig 
abdecken, daS jedes Pixel auf einer solchen Gera- 
den liegt. wobei das Medianfilter 2Mi + 1 Stutz- lO 
stelleo im Abstand a^i hat und dear Pixelwert P des 
mediangefiltertes Bildes MEDIBILD sich eigibt 
sicb aus dem Medianwert Med der Pixelwerte der 
2Mi 1 SttttzsteUen von RIBUJD. 
. c) Bildung des Diflferenzbildes DIFFl s i5 
R1BEI>-MED1BILD fllraUe Np • Np Pixelwerte 
und DurchfOhrung einer Schwellenoperation mil 
der Artefaktschwelle S,rt Bild DIFFl, wobei 
alie Pixelwerte in DIFFl, diebetragsmSBig groBer 
als Salt sind, gleich "fSart Oder "-Sartgesetztwerden, 20 
so dafi das Bild R2BILD entstebt 
d) Bildung des Differenzbildes R3BILD aus dem 
maskierten Ausgangsbild RlBILD und dem Bild 
R2BILD fUr alle Pixelwerte Np • Ni^ 
. e) DurchfUhrung einer zwdten Medianfilterope- 25 
ration in R3BTT J), die voigenommen wird entlang 
einer Vielzahl von durch das Drehzenuum des 
Computertomographen verlaufenden Geraden 
(radiale Linien), welche die BildmaUrix dejardg 
abdecken, dafi jedes Pixel auf einer solchen Gera- 3o 
den liegl, wobei das Medianfilto 2M2+I StOtz- 
stellen im Abstand 8^2 bat und der Pixelwert P des 
mediangefilterten Bildes MED2BILD sich eigibl 
aus dem Medianwert Med der Pixelwerte der ^2 
+ 1 StUtzslellen von R3BILD. 35 
0 Bildung des Dififerenzbildes DIFF2 = 
R1BILD-MED2BIIJD fiir alle Np ■ Np Pixelwerte 
und Durchfiihrung einer Schwellenoperation mil 
der Artefaktschwelle Stn im Bild DIFF2, wobei 
alle Pixelwerte in DIFF2, die beuagsmaBig grfiBer 40 
als San sind, gleich +^'|« oder gesetzt werden, 
so da£ das Bild R4BILD entstebt 

g) Duichftihrung &n& HefpaBtilterung der Pi- 
xelwerte des Bildes R4BILD entlang von Kieis- . 
b6gen urn das Drehzentrum, die als einfache Mil- 4S 
telwertbildung realisiert werden kann, wobei man 
fOr jedes Pixel P im BildR2BILD mit den Kooidi- 
naten xp, yp auf dem durch P verlaufende n Kreis- 
bogen um das Drehzentrum mil Radius yx^y^p 
die Nw + 1 Pixel mit dem Wnkelabstand 50 
-Nw . SPHI. -(Nw+1) • SPm. . . . (Nw- 

1) • SPHI, Nw • SPHI von P sucht und den Wert 
des Pixels P ersetzt durch den Mittelwcrt der 2Nw 
+ 1 Pixel, so daB das Bild WBILD entsteht. 

h) Subuaktion von IBILD und WBILD, wobei 55 
das Eigebnisbild entsteht, in dem die Ringarte- 
fakte beseitigt sind. 

3. Computertomograpfa, der zur Durchfuhnmg des 
Verfahrens nach Anspnich 1 oder 2 ausgebildet isL 
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